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An  approximate  method  In  signal  detection* 
Walter  F.  Freiberger 


Abstract e  A  theorem  from  the  theory  of  Toeplitz 
forms  (ref.  1)  is  applied  to  the  problem  of  esti¬ 
mating  the  best  test  statistic  for  the  detection 
of  Gaussian  signals  in  Gaussian  noise. 


Let  x(t)  be  a  stationary  Gaussian  process,  mean  zero, 
sampled  at  successive  time  points  to  provide  an  observer  with 
a  finite  sample 

2  =  (x^,  x2,...,  xn).  (1) 

The  time  points  are  to  be  closely  spaced  to  produce  almost 
unit  correlation.  It  is  known  to  the  observer  that  (1)  is 
either  a  sample  from  a  Gaussian  ensemble  Qq  with  mean  zero 
and  power-density  spectrum  f0(X)  or  a  sample  from  a  Gaussian 
ensemble  Q-^  with  mean  zero  and  power-density  spectrum  f^CX), 
and  he  is  to  decide  whether  (1)  came  from  £20  or  We  take 

&0  to  be  noise  alone  and  signal  plus  noise  and  write 

f0<X>  =  fn(X)  (2) 

fl<\)  =  fn(x)  +  fg(X)  (3) 

where  f  (\)  and  f  (X)  are  the  power  spectral  densities  of  noise 

XI  5 

and  signal,  respectively}  the  noise  and  signal  processes  are 
here  assumed  independent  and  their  spectial  densities  thus 
additive. 
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Let  H0  denote  the  hypothesis  that  noise  only  is  present 
and  that  xre  observe  signal  plus  noise.  If  Q  =  {u}  :  Qo  U 
denotes  the  sample  space  of  all  possible  realizations  of  the 
process  x(t)  in  a  finite  interval  of  time,  so  that  w  is  a 
function  of  t  <£  (0,T) ,  then  we  define  a  critical  region 


in  the  sense  that 


W  c  Q 


(4) 


if  ,  Hq  is  rejected; 

if  ,  Hq  is  accepted* 

The  probability  that  HQ  be  rejected  though  true  is  denoted  by 
P0(W),  that  it  be  accepted  though  false  by  PjCW*),  and  the 
power  of  the  test,  i.e.  the  probability  of  rejection  of  H0, 
consequently  by 


P1(W)  =  1  -  P1(W*)«  (5) 

It  is  well  known  that  although  perfect  detection  is  not 
possible  in  the  case  we  are  considering,  the  most  powerful 
critical  region  W^p  is  given  by  the  Neyman-Pe arson  test 

%>  =  {x|L(s)  >  c}  (6) 

v/here  the  likelihood  function  L(x)  is  given  by 


L(x)  =  o  x  — rV  • 

povxl,x2,***’xn' 

Here,  p0  (a)  and  p^(oc)  9  the  probability  densities  induced  cn 
Q  by  H0  and  H^,  respectively  are  given  by 


(7) 
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p0<x)  =  7 — ;n72  exP(-  2  (8) 

0  (2it)n/ Vdet  R0  ^  D 

Pi te)  =  7 — exp(-  \  (9) 

(2jc)n/<Vdet  R  -l  ^  1 

R0  is  the  covariance  matrix  of  the  noise,  that  of  the  signal 
plus  noise j  these  are  given  hy 

it 

R0  =  ("2%  J  ^  f0(^-)dXj  V,[i=l,2,...,n}  (10) 

••K 

and  similarly  for  R-^.  Hence 

L(S)  =  K  exp[|  S*(R“1  -  Rj1)*].  (11) 

The  likelihood  function  is  thus  a  monotonically  increasing 
function  of  the  quadratic  form 

^(R"1  -  R”1)^  =  x-Qs  ,  say,*  (12) 

hence 

WMP  =  “  H^1)*  >  c}  .  (13) 

For  practical  purposes  (12)  is  not  a  convenient  expression 
since  its  use  requires  the  inversion  of  large  matrices*  We 
shall  now  show  how  a  theorem  from  the  theory  of  Toeplitz  forms 
(ref,  1)  can  he  applied  here  to  obtain  a  computationally 
feasible  procedure.  A  similar  approximation  has  previously 
(ref.  2)  been  found  effective  for  the  estimation  of  the  spectral 
density  of  a  random  process. 

The  probability  density  of  the  quadratic  form  (12)  can 
be  looked  on  as  being  completely  determined 

(i)  under  hypothesis  H  ,  by  the  eigenvalues  of  the  matrix 

QR0  *  I  -  R£\i  (1W 


If 


(ii)  under  hypothesis  by  the  eigenvalues  of  the 

matrix  . 

QR1  =  -  I.  (15) 


Approximations  to  the  distribution  of  quadratic  forms 
such  as  (12)  are  considered.'  1ft. ref  •  1*  It  is  shown  that  the 
distribution  is  asymptotically  normal,  but  since  the  matrices 
are  nearly  of  Toeplitz  character,  closer  approximations  are 
suggested  by  Toeplitz  theory.  It  can  be  shown  that,  if  p-0 
denotes  the  trace  of  QRq  and  the  trace  of  QR^, 


% 

i  r 

f.,(\)  -  f  A\) 

*s* 

P-o  = 

2%  J 

-X 

X  0 

fjrxr— 

d\, 

(16) 

% 

-i  r 

f,00  -  f„00 

P'1  = 

2%  J 

-L— af- 

dX. 

(17) 

-it 

Hals  suggests  that  we  should  use  for  our  computations  the  quad¬ 
ratic  form 

where 

“  *aJ  eltv"li)Xtrfx7  -  rfx7)a^  Cl8) 

-it  ox 

v,p,  =  1,2, ...,n] 

and  construct  what  we  might  call  the  "approximately  most  povrerful 


test" 

WaMP  =  (sl2^Qa2  >  c}  (19) 


where,  with  (2)  and  (3), 


n  =  afe 


*  n 

;  f  l  xvel 


f  (\) 

fnaufn(0  +  aX- (20) 
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But  we  know  (cf.  ref.  3,  p,  91)  that 

n 


afelJx  *  elvi  ■«». 


(21) 


the  periodogram,  which  is  an  unbiased  and  inconsistent  estimate 
of  the  spectral  density  fOO  of  the  process  x(t).  Thus, 


% 

n 


1  i  f  g  C 

n  X*V  m)  fnix)[fn(0+fs(\)]  dX 


71 


r  fsoo 

=  an  estimate  of  j  f(\)  rTxTTf Jxitf  Jxi'l  dX  *  (22) 

n  n  s 


-7t 


Hence,  the  approximately  best  test  statistic  is  a  weighted 
periodogram.  The  resulting  signal  detection  method  could  there¬ 
fore  be  represented  schematically  as  follows : 


x(t)  — > 


filter 


power  meted  — > 


power 


>c 

<c 


signal  present 
-*>  signal  absent 


The  filter  will  have  to  have  the  relevant  characteristic,  viz 

/ 

^  fnc\)£fn(0+fs()jrj  5 

the  power  meter  is  a  quadratic  integrator.  Usually,  shapes  of 
the  spectral  densities  will  be  somewhat  as  follows: 
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With  this  procedure,  then,  ve  do  not  have  to  invert  large 
matrices  but  can  construct  simple  physical  devices j  it  is  also 
suitable  for  high-speed  computation. 

To  test  the  accuracy  of  this  method,  the  frequency 
functions  corresponding  to  the  most  powerful  test  based  on 

x*Qx 

and  the  approximately  most  powerful  test  based  on 

x*Qax 

were  computed  and  compared  numerically  for  a  simple  case. 

We  chose  for  the  covariance  matrix  of  the  noise  the  unit 

matrix 

R0  =  N  =  I 

and  for  the  covariance  matrix  of  the  signal 

Rl  -  Rq  =  S  =  v,p.  =  l,2,,.,,n} 

with 


SV(i 


e^-^ax  =  pl^l 

ll-pe^l2 


The  spectral  density  of  the  noise  is,  hence,  fn(\)  =  1  and 
that  of  the  signal 


fsc\>  = 

is  the  Poisson  kernel  which, 


— J. gfil .... 

Il-pe^j2 

for  values  of  p  fairly  close  to  1, 


is  narrow-band,  ^he  parameters  chosen  were 


n  —  20  ,  p  —  o7  • 

If  we  denote  the  eigenvalues  of  S  by  Xy,  those  of 
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Q  =  N"1  -  (N+S)"1  =  I  -  (I+S)"1 


will  be 

Ue  denote  by  gQ(x) 
eigenvalues 


the  frequency  function  corresponding  to  the 

1  -  ifc; 


and  by  g-^Cx)  that  corresponding  to  the  eigenvalues  Xy,  These 
are  then  the  ’’exact"  frequency  functions* 

For  the  approximate  theory,  we  define  the  matrix  Q  with 

cl 

elements 


% 


-  msLm 

2n 


i(v-|x)  X  . 

fn(n)lf“(\)+fs(X)j  e  dX 


-K 


a  JL 

2it 


it 


1-  9C 


Ps1Kl2 

1  1  +  1  -  K 

j1  H-pe^l2 

~x 


i  jl  _ 

2lt  Jl-pe1^!2  +  1 


-s  dX  . 
P 


-7t 

By  choosing  the  constant  r  as  the  solution  of  the  equation 
r[(l-p2)  +  (1+p2)]  =  p(l+r^) 
such  that  0  <  r  <  can  be  brought  to  the  form 


The  eigen- 


•'vn  2%p 
-it 


ei(v  ^)XdX 


1 1-re 


iXi2 


it 


p(l-p2)  2% 

=  _ E. . 

p(l-p2) 


1=3:  ei(v“^Xd\ 


ll-re^  |2 


•  r(v-^  , 


since  r 


( v—  Jl) 


is  the  Fourier  transform  of 


JdL 


values  of 


1-re 


iX72 


Qa  =  { ^vp.  5  vjti  =  !>•••}  n}j  n-20j  P  -  »7 

are  computed  to  give  the  frequency  function  Y0(x)  which  is  to  he 
an  approximation  to  gQ(x)  as  defined  above* 

Similarly,  the  eigenvalues  of 

«aRl  =  VRo  +  S)  =  V1  +  » 

are  found  and  used  to  determine  the  frequency  function  Y^Cx) 
which  serves  as  an  approximation  to  g-^(x)  as  defined  above. 

Figures  1.  and  2  show  the  curves  gQ(x),  Y0(x)  and  gj(x), 
Y^x),  and  agreement  is  seen  to  be  very  good. 


Acknowledgment .  Thanks  are  due  to  Professor  Ulf  Grenander  on 
whose  work  and  suggestions  this  paper  is  based. 


References 

it  _ 

1.  U.  Grenander  and  G.  Szego,  "Toeplitz  Forms  and  their 

Applications",  University  of  California  Press,  1958. 

2.  W.  Freiberger  and  U.  Grenander,  Approximate  distributions 

of  noise  power  measurements,  Quart.  Appl*  Math.  17, 
272-283  (1959) • 

3.  U.  Grenander  and  11.  Rosenblatt,  "Statistical  Analysis  of 

Stationary  Time  Series",  John  Uiley  &  Sons,  New  York$ 
Almquist  &  Uiksell,  Stockholm,  1957* 


9 


9 

9 


Csl 

fO 


oo 


o 


FIGURE 


